Cluster Dual Fermion Approach to Nonlocal Correlations 
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We formulate a general cluster Dual Fermion Approach to nonlocal correlations in crystals. The 
scheme allows the treatment of long-range correlations beyond cluster DMFT and nonlocal effects in 
realistic calculations of multiorbital systems. We show that the the simplest approximation exactly 
corresponds to free cluster DMFT. We further consider the relation between the two-particle Green 
functions in real and dual variables. We apply this approach by calculating the Green function of 
the Hubbard model in one dimension starting from the two-site cluster DMFT solution. The result 
agrees well with the Green function obtained from a DMRG calculation. 

PACS numbers: 



(N 



I 

o 
o 



> 

(N 
O 

o 
o 



I. INTRODUCTION 

One of the successful routes to the description of 
strongly correlated systems is related to the Dynami- 
cal Mean Field Theory (DMFT>i.. In this scheme the 
system is mapped onto an effective local quantum im- 
purity problem in a self-consistently determined bath. 
The self-energy in the DMFT approach is local in space 
but frequency dependent. However, there are many phe- 
nomena for which non-local correlations are important 
and often correlations arc long-ranged. Among these are 
Luttinger-Liquid formation in low-dimensional systems, 
non-Fermi-Liquid behavior due to van-Hove singularities 
in two dimensions or d-wave pairing in high- Tc super- 
conductors. Obviously, DMFT is not suitable for the 
description of such systems and there are cases in which 
DMFT even fails qualitatively, as is for example the case 
for the one-dimensional Hubbard model. For the descrip- 
tion of these non-local correlations it is desirable to com- 
bine short-range cluster many-body physics, like forma- 
tion of singlets, and long-range correlations. Recently the 
novel Dual Fermion approach for the treatment of nonlo- 
cal correlations has been developed^ . Here we formulate 
a general cluster (or multi-orbital) scheme for non-local 
correlations. 

Frequently used approaches to account for non- 
local correlations beyond DMFT comprise the clus- 
ter approaches or the so-called Dynamical Cluster ap- 
proximation (DCA) in k-space^, real space periodio^ 
and free cluster approaches^, as well as the Cellular- 
DMFT6(CDMFT). In the latter approach, the single-site 
impurity of the DMFT is replaced by a cluster of impu- 
rities. The CDMFT scheme however is restricted to rela- 
tively small cluster sizes due to computational feasibility 
and only short-range correlations can be treated within 
this approach^. 

Recently, steps have been taken to go beyond DMFT 
and to treat long-range correlations. One of them 
is the Dynamical Vertex Approximations and similar 
approaches^iiS, where a diagrammatic expansion around 
DMFT is made. A principally new scheme with a fully 



renormalized expansion called Dual Fermion Approach 
has been proposed-S, which is based on the introduction 
of new variables in the path integral representation. Sim- 
ilar schemes for the strong coupling expansion have been 
discussed in terms of Hubbard operatorsii. 

Here we present a cluster generalization of the Dual 
Fermion approach. This allows the treatment of clusters 
or multiorbital atoms within the Dual Fermion frame- 
work and can describe long-range correlations in realistic 
systems. The paper is organized as follows: first, we 
present the formalism and then show that for the case 
of non-interacting dual fermions one exactly obtains the 
self-consistency equation of the CDMFT. In the next sec- 
tion we consider the relationship between the two-particle 
Green functions in real and dual variables. Then we de- 
scribe the calculation procedure. We finally apply our 
approach to the Hubbard model in one dimension start- 
ing from a two-site CDMFT solution. 



II. FORMALISM 

Our goal is to find an (approximate) solution to the 
cluster lattice problem described by the imaginary time 
action which in cluster or multiband notation reads 



S[c*,c] 



(1) 



Here i?o-(k) is the one-electron part of the Hamiltonian, 
u! — {2n + l)7r//3,n = 0,±1,... are the Matsubara fre- 
quencies, /3 and fi are the inverse temperature and chem- 
ical potential, respectively, a =f , i labels the spin projec- 
tion and c*, c are Grassmannian variables. The indices i 
label the translations of the super-cell and the k- vectors 
span the reciprocal lattice in the reduced super-cell Bril- 
louin zone. Here we take into account a general type of 
interaction, Hi^f It is important to note that it can be 
any type of interaction inside the multiorbital atom or 
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cluster. The only requirement is that iJjnt is local. For 
example, the general Coulomb interaction has the form 



where 



Hint[c*,Ci\ 



-y 



dr Ui234clclc4C3 , 



(2) 



where U is the general symmetrized Coulomb vertex and 
e.g. 1 = {wimiCTi} comprehends frequency-, orbital- and 
spin degrees of freedom and summation over these states 
is implied. In order to capture the local physics we in- 
troduce a cluster impurity problem just in the spirit of 
CDMFT in the form 



imp [ 



, c = - > c, 



-I- Hint[c*,c] , 



(3) 



where A is an as yet unspecified hybridization function 
describing the interaction of the impurity with an elec- 
tronic bath. We suppose that all properties of the im- 
purity problem are in principle known, i.e. the single- 
particle Green function ga-u) and the irreducible vertices 
7^'*'', 7^, etc. are known. Our goal is to express the Green 
function G^^k and vertices 7 of the original lattice prob- 
lem via these quantities. 

Since A is local, one may formally rewrite the original 
lattice problem in the following form: 



■> imp ['-■i^icr J ^wiCTj 



5[c*,c] = Sn 
i 



(4) 



We introduce spinors c^ka = (• ■ • , c^kam, ■ ■ ■), c*^^ = 
(. . . , c*j.g.^, . . .). Omitting indices, in matrix- vector no- 
tation, the Gaussian identity that facilitates the trans- 
formation to the dual variables is 



J exp (-rif - rBc - c*Bf) V[i, f * 
det(i) exp(c*BA-^Bc 



(5) 



which is valid for arbitrary complex matrices A and B. 
In order to decouple the non-local term in Eqn. 31 we 
choose 



A = .gji(A^,-iJ,(k))-igji , 

B = 9Zl , (6) 



where g^ja is the Green function matrix in the orbital 
space (m,m') of the local impurity problem. Using this 
identity, the lattice action can be rewritten in the form 

5'[c*, C, f*, f] = ^ ' 'S'sitc,;: + 
i 

(A^.-i/.(k))-^ 5jjLk.] , (7) 



i i 

(8) 

Here the summation in the last term over states labeled 
by k has been replaced by the equivalent summation over 
all sites. The Gaussian identity can further be used to es- 
tablish an exact relation between the lattice Green func- 
tion and the dual Green function. To this end, the parti- 
tion function of the lattice is written in the two equivalent 
forms 

Z = y cxp(-5[c*,c])2?[c,c*] = 

exp(-5[c*,c,r,f])P[f,r]P[c,c*] , 



where 



Zf - n "^^^ (A^^ - i?a(k)) g^„] 



(9) 



(10) 



By taking the functional derivative of the partition func- 
tion, Eqn. [51 with respect to the Hamiltonian, i.e. 



1 



5Z 



(11) 



one can obtain the following exact relationship between 
the dual and lattice Green functions: 

G^ka - (5a.a(A^a-i?.(k)))"'G5ka X 

X ((A^,-i7,(k)) V(A„a-i/.(k))~^ , 

(12) 

where the lattice Green function is defined via the imag- 
inary time path integral as 



Gi2--| j cic^exp(-5[c*,c])I?[c,c*] 



(13) 



and similarly for the local Green function g and dual 
Green function G'' with Z and S replaced by the corre- 
sponding expressions. 

We now wish to derive an action depending on the dual 
variables only. This can be achieved by integrating out 
the original variables c,c*. The crucial point is that this 
can be done for each site separately: 

j exp(-S'sito[c*,c,,f*,fi])X'[ci,c*] = 

This equation can be seen as the defining equation for the 
dual potential y[f*, f]. Since S^itc contains the impurity 
action, expanding the remaining part of the exponential 
and integrating out the original variables corresponds to 
averaging over the impurity degrees of freedom. Equating 
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FIG. 1: The first two lowest order diagrams for tlie dual self 
energy E**. 



the resulting expressions by order, one finds that the dual 
potential in lowest order approximation is given by 



(15) 



where 



7l234 — 5ll'522' 



imp _ (O)imp 
Al'2'3'4' Al'2'3'4' 



33'3ff4'4 ' 



(O)imp 

Xl234 = 514523 - 513524 



(16) 



is the fully antisymmetric irreducible vertex and the 
local two-particle Green function of the impurity model 
is defined as 



imp 

Xl234 — 



1 



-'imp 



C1C2C3C4 exp (-S'inipic*, c]) V[c, c*] 



(17) 

The dual action now depends on dual variables only 
and can be written as 



jkcr 



The bare dual Green function is given by 



G 



k(T ^ 5WCT 



(18) 



(19) 



For the nonlocal part of the self energy defined by 
Snonioc = S — Sloe wc find a simple matrix relation to 
the dual self-energy: 



^— 1 y]~' 

^nonloc d 
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(20) 



The local part of the self energy is obtained by the so- 
lution of the impurity problem. In order to obtain the 
nonlocal contribution, we thus need to calculate the dual 
self-energy. This is achieved by performing a regular di- 
agrammatic series expansion of the dual action, Eqn. 
[T51 and considering the lowest order diagrams for Ed, 
constructed from the irreducible vertices and the dual 
Green function as lines. The diagrams considered here 
are shown in Fig. [TJ The lowest order diagram is local 
while the next diagram already gives a nonlocal contri- 
bution to the self energy. 

So far we have not established a conditon for A. We 
require that the first diagram in the expansion of the 
dual self-energy should be equal to zero at all frequencies. 
Since 7''*' is local, we can use the following condition: 



E 



Gd 







(21) 



In the simplest approximation, which corresponds to non- 
interacting dual fermions, the full dual Green function is 
replaced by the corresponding bare Green function and 
the above condition can be reduced to 



E 

k 







(22) 



which is equivalent to the self-consistency condition for 
the hybridization function in free-cluster CDMFT— i^i^. 
The self-consistency condition for A which satisfies this 
condition is 



Did ■ 



5 ^G'focG'loJ 



(23) 



If we restrict the matrix A^^o- to be equivalent for all 
cluster sites with periodic boundary conditions then this 
approximation leads to the DCA scheme^. 



III. TWO-PARTICLE EXCITATIONS 

In order to find the exact relation between the four- 
point correlation function in dual and conventional vari- 
ables we have to calculate the second derivative of Z 
with respect to iJcr(k)^A, Ha{k.)pi, using the two equiv- 
alent representations of the partition function, Eqn. [Sj 
Differentiating the partition function in its original form 
gives us 
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52 Z 



Z 5H^x5Hp, 



(24) 



while differentiating the second expression for Z using 
the above relation for Green's function, Eqn. [121 yields 
for x>^tJ.^p after some straightforward algebra: 

[(A - ® [(A - H)-^g-^Gag-HA - H)-^\\ + 
[[(A - li)-^g-^Gag-\A - H)-'] ® (A - H)-'] ^^^^ + 
[(A - H)-'g-'] [(A - Hr'g-'] ^„ xLnr x 
x[g-\A-H)-']^^^[g-\A-H)-']^^. (25) 

Here xL«r = (^fifLfnf}) is the dual four-point corre- 
lation function and 

{A B)\pi,p = A\pByp — A\pB^p 

is the antisymmetrized direct product of two matrices. 

As one can see from the two-particle excitations 
for dual fermions coincide with those for real fermions. 
Thus to get the information about the instabilities in the 
system under consideration it is sufficient to sum up the 
ladder for the two-particle dual fermion Green function 
in the ladder approximation: 

d _ do I do d (OPi) 
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IS 



DMFT 



outer loop 



g.7 



(4) 



inner 
loop 



— ' 



executed until self-consistency. Self-consistency for both 
loops is usually reached after a few iterations (depending 
on the system). The computational cost for the calcula- 
tions aside from DMFT is less than for the DMFT itself, 
whereby the computation of the vertex is the computa- 
tionally most expensive part. We use the continuous-time 
quantum Monte Carlo impurity solver^^ for the solution 
of the impurity problem and for the calculation of the 
irreducible vertex. 



V. APPLICATION TO ID HUBBARD MODEL 

We consider the one-dimensional Hubbard model de- 
scribed by the following Hamiltonian: 



Ho + Hint = t ^ cl_^_^Ci + U'^ n-i^Tiii 



(28) 



FIG. 2: Calculation procedure 



or 



d _ do I do d (9,7) 

where the first equation is written for the particle-particle 
channel and the second for the particle-hole one, and 
^do — jg bare two-particle dual Green func- 

tion. 



For the half-filled case the main physics is related with 
the formation of a local singlet pair. Therefore we start 
from a two-site cluster DMFT solution and then include 
the long-ranged non-local effects via the Cluster Dual 
Fermion approach. When this one-dimensional system is 
treated as a chain of two-site clusters as depicted in Fig. 
[3l the tight-binding Hamiltonian for this model is readily 
shown to be 



Ho{k) 







t{l 



-ika\ 



(29) 



IV. CALCULATION PROCEDURE 

The calculation procedure is as follows: Starting from 
a starting guess for A, e.g. A = 0, we solve the im- 
purity problem and obtain a new local Green function 
g. From this we calculate the bare dual k-dependent 
Green function via Eqn. In order to reach the self- 
consistent free-cluster DMFT solution we calculate the 
local part of the lattice Green function via Eqn. [T^] and 
insert it into Eqn. [53] to obtain a new A with which we 
again solve the impurity problem. This loop, which is 
closed by the dashed arrow in Fig. O is repeated un- 
til the self-consistent CDMFT solution is reached. With 
the thus obtained A we then calculate the irreducible 
vertex 7^"'^ Now we do not follow the path indicated 
by the dashed line but calculate an approximation to the 
dual self-energy by summing up the first diagram(s) in 
the perturbation series expansion. From this and the 
bare dual Green function an approximation for the dual 
Green function GJ^j^.^ is obtained which is subsequently 
used in the diagrams. This inner loop is executed until 
self-consistency. The dual Green function and the lat- 
tice Green function are then again used to obtain a new 
hybridization function, which serves as input for the cal- 
culation of a new local Green function and renormalized 
vertex in the impurity solver step. This outer loop is also 
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FIG. 3: Schematic representation of the ID chain as a chain 
of two-site clusters. 

Due to the absence of a Mott transition in one dimen- 
sion, the system is an insulator for any finite value of the 
on-site repulsion U. It is known that single-site DMFT 
gives a qualitatively wrong answer, i.e. predicts the sys- 
tem to be metallic even for U as large as J7/t = 6. Here we 
compare our results with the one obtained by a Density 
Matrix Renormalization Group (DMRG) calculation^-, 
since DMRG is known to reproduce the spectral prop- 
erties of ID systems quite well. The free cluster DMFT 
solution obtained in our calculation is qualitatively cor- 
rect and reproduces the DMRG solution quite well, since 
short range fluctuations between nearest neighbours are 
predominant. For our calculations we use the parameters 
U/t = 6 and /3 = 20. The DMRG solution corresponds 
to T = 0. For the calculation we considered only the first 
two lowest order diagrams in the series expansion shown 
in Fig. [TJ The result of our calculation is depicted in Fig. 
m where we show the imaginary part of the on-site Green 
function as a function of Matsubara frequencies on one of 
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FIG. 4: (color online) Local Matsubara Green function on 
the Matsubara axis obtained from DMRG for T = in com- 
parison with the results obtained from DMFT and from fully 
self-consistent dual fermion calculations. For the 2-site free 
cluster DMFT the renormalization of the vertex has a small 
effect since the CDMFT is already a good starting point. For 
the single site calculation the renormalization is essential since 
DMFT even gives a qualitatively wrong answer, while the dual 
fermion result correctly predicts the system to be an insulator. 



out changing the vertex. We further plot the result for 
a single-site calculation. Here single-site DMFT wrongly 
predicts the system to be metallic. Hence the DMFT 
solution strongly differs from the exact solution and the 
renormalization of the vertex becomes important. Never- 
theless, the dual fermion solution converges to a metallic 
solution within a few iterations. 



VI. CONCLUSIONS 

We have generalized the recently proposed Dual 
Fermion Approach to the multiorbital case, facilitating 
the treatment of multiorbital systems within this frame- 
work. We further established the relation between the 
four-point correlation functions in real and dual vari- 
ables and found that the two-particle excitations for real 
fermions and dual fermions coincide. The approach was 
applied to the one-dimensional Hubbard model starting 
from the free two-site cluster DMFT solution. Although 
the CDMFT solution already quite well reproduces the 
DMRG solution, the cluster dual fermion solution con- 
siderably improved this result. The cluster formulation 
allows to combine this approach with realistic density 
functional calculations and thus opens a new way to de- 
scribe long-range correlations in real systems. 



the cluster sites. In this calculation the renormalization 
of the vertex has a small effect, since the CDMFT result 
is already close to the exact solution. However, the free 
cluster DMFT result is considerably improved and the re- 
sult agrees very well with the one obtained in the DMRG 
calculation. Almost the same result was obtained with- 
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